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Investigation of the spacelike and timelike electromagnetic form factors of hadrons, 
within a relativistic microscopical model characterized by a small set of hypothesis, could 
shed light on the components of hadron states beyond the valence one. Our relativistic 
approach has been successfully applied first to the pion and then the extension to the 
nucleon has been undertaken. The pion case is shortly reviewed as an illustrative example 
for introducing the main ingredients of our approach, and preliminary results for the 
nucleon in the spacelike range —10 (GeV/c) 2 < q 2 < are evaluated. 

1. INTRODUCTION 

The global analysis of the electromagnetic (EM) form factors in both the spacelike (SL) 
region and the timelike (TL) one, within an approach based on the Light-Front language 
[EH I2j, allows one to investigate in great detail the issue of non- valence components of 
both hadron and photon wave functions. In Refs. [HU H], the pion EM form factor, in 
the range of momentum transfer —10 (GeV/c) 2 < q 2 < 10 (GeV/c) 2 , has been analyzed 
within our relativistic approach, based on: i) a microscopic Vector Meson Dominance 
(VMD) model for the dressed vertex of the photon, and ii) a simple parametrization for 
the emission/absorption of a pion by a quark. Indeed, the pion represents a suitable 
introduction for illustrating the main features of our approach [ OH E] . 

Light-front dynamics (LFD) approaches [ |2] yield a unique possibility to study the 
hadronic state, in both the valence and the nonvalence sectors, since within LFD a Fock 
series for the state of a fermionic system, like 



\meson) = \qq) + \qqqq) + \qq g) 

\baryon) = \qqq) + \qqq qq) + \qqq g). 
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is a meaningful tool for investigating hadronic properties. It is worth noting that the 
so-called zero modes allow one to produce the breaking of the chiral symmetry through 
a different mechanism from the standard one, so that the LF vacuum acquires a more 
simple structure. Moreover, coupling the study of the EM form factors in the TL region 
with the analysis in the SL one, yields the possibility to address the vast phenomenology 
of hadronic resonances, e.g., the vector meson (VM) propagation, and consequently to 
impose strong constraints on dynamical models devised for a microscopical description of 
hadrons. Finally, one can obtain insight into the two-body currents associated to the qq 
pair production, that plays a very relevant role in a reference frame where q + 7^ (see, 
e -g- [El)) namely the frame mandatory for carrying on our global analysis of the EM form 
factors. 



2. The Mandelstam formula for the EM current 

Our guidance for constructing a microscopical model that embeds, at large extent, the 
advantages of a LF approach is given by the Mandelstam formula [E] for the matrix ele- 
ments of the EM current of a system of interacting constituents. In a global investigation 
of SL and TL regions one needs to change frame, from the q + = frame (the standard 
choice within LFD, that always leads to q 2 < 0) to a q + ^ frame, that allows both 
positive and negative q 2 . Therefore, a covariant expression for the matrix elements of the 
EM current of hadrons is a useful tool for constructing viable approximations. 

For instance, in the TL region the Mandelstam formula reads as follows 

f = -«e / [S Q (k Q )\ h (q - k, k Q , P h )S(q - k)T»(k, q)S(-k)X h (k, k Q , P~ h )} (1) 

where S(p) = 1/(j4 — m + te), with m the mass of the constituent struck by the virtual 
photon, Sq(p) is the propagator of the spectator constituent (e.g. a quark, or a diquark 
in a simple picture of baryons), Xh{k, k', Ph) is the hadron vertex function, which contains 
the Dirac structure, i.e. a proper combination of Dirac matrices, and a momentum de- 
pendence; P£ {Pj*) is the hadron (antihadron) momentum; T^ 1 {k ) q) is the quark-photon 
vertex, and g M is the virtual photon momentum. The corresponding SL expression can 
be obtained by applying the suitable changes to the hadron momenta and to the Dirac 
structure. The pion case is illustrated in detail in Refs. [El El- It should be pointed out 
that the momentum components of both and T M are necessary for regularizing Eq. 

Projecting out the Mandelstam formula on the Light-Front hyperplane, through a k~ 
integration, allows us to introduce sensible approximations for constructing our phe- 
nomenological model. Let us remind that, presently, solving the Bethe-Salpeter equation 
for fermionic system represents a big challenge, and therefore phenomenological studies 
could play a very useful role. 

Our assumptions for evaluating the k~ integration of Eq. are: i) the momentum 
components of the vertex functions, both for the pion and for the vector mesons entering 
in our VMD model for r M (/c,g), vanish in the k~ complex plane for \k~\ — > 00, ii) the 
contributions of the vertex function singularities can be neglected. These approximations 
lead to the following issues: i) how to connect the Fock language with the Bethe-Salpeter 
one, e.g. how to describe both the amplitude for the emission/absorption of a pion by a 
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quark (i.e. the nonvalence component of the pion wave function) and the gg-pion vertex 
(i.e. the valence component), ii) how to model the quark-photon vertex. In the limit 
m v — > 0, only the contribution generated by the nonvalence component of the hadron 
state, i.e. a higher Fock component, is acting in both TL and SL regions [ |3J Hj, and 
therefore the analysis is remarkably simplified. In particular, in order to describe the 
emission/absorption of a pion by a quark, following Ref. [ 7] we assumed a pointlike 
interaction. It is worth noting that in the present approach, the coupling constant is fixed 
by the charge normalization of the pion form factor. 

As to the quark-photon vertex, the bare term can be dropped, given our assumption on 
the pion mass, and one remains with the effect due to the qq production. In view of this, 
a VMD approximation has been applied, namely the relevant plus component of r M (fc, q) 
is spanned over a VM basis as follows 



r+(*, q) = V2j2 [ex- UK k - P n )\ K{k, P n ) - X X l f Mv i *v\ (2) 

n,A L J (<? 2 - Ml + iM n T n {q 2 )) 



where fy n is the decay constant of the n-th vector meson into a virtual photon (to be 
calculated within our model), M n the corresponding mass, T n (q 2 ) the total decay width 
(equal to T n q 2 /M 2 for q 2 > and vanishing for q 2 < 0) and e\ the VM polarization. The 
VM vertex amplitude is e A ■ V n (k, k — P n ) A n (k,P n ), with a proper Dirac structure, 

V n (k, k — P n ), that generates the Melosh rotations for 3 S*i states and a momentum- 
dependent part indicated by A n {k,P n ). In the kinematical region corresponding to the 
valence sector, the Bethe-Salpeter amplitudes, for both pion and vector mesons, are esti- 
mated through the eigenvectors of a relativistic simple model for pseudoscalar and vector 
mesons [E]. 

An estimate of the probability of the valence component, P q q, is necessary for a proper 
normalization of such a component for both pion and vector mesons. In Ref. [0] a simple 
and intuitive model for P qq is presented. 

In the actual calculation, 20 vector mesons are taken into account to reach convergence 
up to \q 2 \ = 10 (GeV/c) 2 . The corresponding eigenstates are calculated from the model 
of Ref. [ 9j, while for the masses a different choice has been made: for the first four 
vector mesons, the experimental values (together with the corresponding widths) have 
been considered, while for M n > 2.150 GeV the model eigenvalues have been adopted. 
We have introduced two adjusted parameters: i) a single width, T n , for the vector mesons 
with mass M n > 2.150 GeV, ii) the relative weight, wvm, of the so-called "instantaneous" 
contributions of the pion and the vector mesons. These instantaneous terms are generated 
by the instantaneous part of the quark propagator, and they play the central role in the 
evaluation of the pion form factor, in the limit of a vanishing pion mass. We have chosen 
for the experimentally unknown widths a value T n = 0.15 GeV, similar to the experimental 
width of the first four VM's; while for wvm we have found that wvm = —0.7 is a good 
value for a global fit (see Fig. 1), and Wvm = —1-5 gives an improved description of the 
p peak region (see Figs. 10 and 11 of [ 3 ). 
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Figure 1. Pion EM form fac- 
tor in the SL and TL regions 
vs q 2 . Solid line: calculations 
with the pion wave function 
from the model of Ref. [Ej, 
and a relative weight wvm = 
—0.7 for the "instantaneous" 
terms. Dashed line: the same 
as the solid line, but with the 
asymptotic pion wave func- 
tion [EH- Data from the col- 
lection of Ref. [EDI. (After 
Ref. [El) 



3. The EM nucleon form factors 



The evaluation of the Sachs form factors of the nucleon largely proceeds along the 
same line of the pion. We consider the nucleon as a system of quarks with mass equal 
to .220 GeV, and then we model the quark-nucleon amplitude in analogy with the pion 
case. First of all, we adopt a Dirac structure of the quark-nucleon vertex as suggested by 
an effective Lagrangian density like the one of Ref. [Ill], viz 

C eff (x,T 1 ,T 2 ,T 3 ,T N ) = e abc J d 4 x 1 d A x 2 d 4 x 3 J r (x 1 , x 2 , x 3 , x) x 



a 



x 3 ) 



)V 3 T (-. &>) 



n,T2,T3 



i) N (x)T TN + 



+cyclic + h.c. 



(3) 



where T{x\ y scalar distribution that properly weights the quark coordinates; 

e a b c takes care of the color degrees of freedom; is the mass of the nucleon; q(x,r) = 
q(x) T T is the field that creates a quark u or d depending upon r; the charge conjugated 



field is q c (x,r) = Cq(x,r)& = Cf{x) T T C = Cq 1 



ix) i T„ 



T* = i'j 2 'j q T (x)i T y T* and 



has a parity opposite to the quark field (this entails a pseudoscalar or an axial operator 
between q a and q b c , in order to recover the correct intrinsic parity of the nucleon). The 
isospin field = i r y T* creates a charge conjugate isospin state, namely i r y x* 2 - For 
the present time a = 1, i.e. we have considered a pseudoscalar coupling only. 

In what follows, given the preliminary character of this presentation, we deal with the 
discussion of the EM form factor in the SL region. The matrix elements of the macroscopic 
EM current of the nucleon, 



(a',p'\r \p,a) = U N (p',a') 



-F 2 (Q 2 



p + p M 
2M N 



+ (F 1 (Q 2 ) + F 2 (Q : 



r 



U N (p,a) (4) 
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where Q 2 = —q 2 and Fu2) is the Dirac (Pauli) nucleoli form factor, are approximated 
microscopically as follows 

{a',p'\f \p,a) = N c / — — / — - x 



.. (27r) 4 7 (2tt) 4 

Tr^ihMX^') S-\h) S-\k 2 )T% ^(k^k^k^p)} (5) 

where the factor N c is a color weight, the trace is performed over isospin and Dirac ma- 
trices, $jy-(fci, k 2 , ks,p) is the quark-nucleon Bethe-Salpeter amplitude [H2], that contains 
both a Dirac structure as suggested by Eq. (j3J) and a momentum component, as in 
the pion case; X M is the quark-photon vertex with isoscalar and isovector contributions, 
namely 

/ 1 _i_ T \ /i t \ -l- T M tm _ T M 

Z* = (^U^+fV^U^^^+^^T^-^^ + r^ (6) 



2 / " V 2 / a 2 2 

The terms and Ij V contain a pointlike (bare) contribution and a VMD contribution, 
with isoscalar (a>-like) and isovector (p-like) content, respectively, viz 

^ siIV) (k,q)=Af IS{ iv)0(p + -k + )6(k + )r + 6(p' + - k + ) 6(k+ - p + ) x 

{z b Mis^-f + z v Klk, q, is (iv)]} (7) 

where A/75 = 1/6 and A/rs = 1/2, the quantities Z& (bare term) and Zy (VMD term) are 
unknown renormalization constants to be extracted from the phenomenological analysis 
of the data. The theta functions, acting on the plus component of the momenta involved 
in the process, implement the kinematical constraints obtained from the projection of 
Eq. (j3J) on the LF hyperplane, as in the case of the pion. In particular, they identify 
the kinematical region relevant to the physical processes depicted in Fig. 2, that can be 
interpreted as valence and nonvalence contributions. Moreover, the vertical dashed lines 
in Fig. 2 allow one to count how many constituents are in flight, at a given LF time. 
It is worth noting that while the VMD is uniquely associated to the pair production by 
the virtual photon, (cf contribution (b) in Fig. 2), i.e. to the nonvalence component of 
the nucleon, the bare term (7^) can yield a contribution both in the elastic quark-photon 
channel (or triangle diagram, cf contribution (a) in Fig. 2) and in the pair production 
process. 

The momentum dependence of the quark-nucleon amplitude in the valence sector, where 
the spectator quarks are on their-own A; _ -shell, has been modeled adopting a power-law 
Ansatz [ 13 , like 

W N (ki, h, k 3 ,p) ~ 1 o (8) 

[^ + M 2 (l,2,3)f 

where (3 has been fixed through the anomalous magnetic moments of the nucleon. In 
particular, with (3 = 0.118 GeV, one obtains fi p = 2.878 (Exp. 2.793) and \i n = —1.859 
(Exp. —1.913), and Mo(l,2,3) is the LF free mass for a three-quark system. In the 
non-valence sector, relevant for evaluating the Z-diagram contribution, the momentum 
dependence of the quark-nucleon amplitude is approximated by 

GN(klMk3 ' P) ^[P 2 + mi,2)r {^ + Mo 2 (3',2)] + [^ + 4(3M)]} (9) 
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Figure 2. Diagrams contributing to the nucleon EM form factors in the spacelike region. 
Diagram (a) (triangle or elastic term) illustrates the valence sector contribution, where < 
kf < P^ (i = 1, 2, 3), and the symbol x on a quark line indicates that the corresponding 
quark is on its fc~-shell. Diagram (6) shows the nonvalence sector contribution, P£ < 
k£ < P'%, due to the pair production. 



where Mo(i,j) is the LF free mass for a two-quark system. This Ansatz is different from 
the one used for the pion case, where a simple constant vertex was adopted, since the 
sixth-fold integration needs a stronger regularization. 

For the present, preliminary calculations in the SL region we have the following adjusted 
parameters: i) Zf, and Zy, that are renormalization constants for the pair production 
terms in the quark-photon vertex (cf Eq. (J7J)), ii) two parameters for a running mass in 
the valence sector (with a range of variation of about 10%, for —10 GeV 2 < q 2 < 0), iii) 
two parameters for a running mass in the non- valence sector (with a range of variation of 
about 20%, for -10 GeV 2 < q 2 < 0). 

The Sachs form factors can be obtained from the macroscopic current by applying the 
proper traces to Eq. (|3J) and then by using the microscopical evaluation of the relevant 
matrix elements. In particular in the frame where p' ± = p± = and q + = y/Q 2 

Q 2 



4M 2 N 



M 2 N 
2p+p H 



-Tr 



N 



2M 



N 



F 2 N {Q 



2 y +P H 



2M 



N 



F?(Qr) + F»{Q 
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N 



2M, 



-7 



TV 



r M 



(Q 



Q 2 



F 1 (Q 2 )+F 2 (Q 2 
f f + M N 



2M 



N 



~ F 2 (Q 



/I i 1 
,2^P +P 



2M 



N 



+f 2 n (q 2 : 



1 



j) + M N / 
2M, 



(10) 



N 



In Figs. 3 and 4, the ratio G p E (Q 2 )fx p /G p M (Q 2 ), and the form factors G P M (Q 2 ), G n E {Q 2 ), 
G^jlQ 2 ) are shown. One of the main feature of our preliminary calculation is the ex- 
planation of the possible zero in the ratio G p E {Q 2 )ii p /G p M {Q 2 ) in terms of the interplay 
between the Z-diagram contribution, i.e. higher Fock components, and the elastic contri- 
bution. Such an interpretation appears supported by the encouraging agreement between 
theoretical calculations and experimental data. 
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Figure 3. Left panel. The ratio G p E (Q 2 )fi p /G p M (Q 2 ) vs Q 2 . Solid line: full calculation, 
corresponding to the sum of all the contributions to G P E {Q 2 ) and G P M (Q 2 ), i.e. triangle plus 
pair production terms. Dotted line: triangle elastic contribution to G P E {Q 2 ) and G P M {Q 2 ). 
Data: from [Ej. Right panel. G P M {Q 2 ) / \i v G d(Q 2 ) vs Q 2 , with the same notations as in 
the left panel, and G D (Q 2 ) = 1/(1 + Q 2 /0.71) 2 . 




Figure 4. Left panel. G E (Q 2 ) vs Q 2 . Solid line: full calculation, corresponding to 
the sum of all the contributions to G E {Q 2 ), i.e. triangle plus pair production terms. 
Dotted line: triangle elastic contribution to G E {Q 2 ). Data: from [E]. Right panel. 
Gm(Q 2 ) I l^nGniQ 2 ) vs Q 2 , with the same notations as in the left panel, and Gd{Q 2 ) = 
1/(1 + Q 2 /0.71) 2 . 
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4. Conclusions & Perspectives 

In [ E] , we constructed a microscopical model for the pion form factor, within a LF 
approach, obtaining a very nice description of experimental data. The development of 
this model for calculating EM nucleon form factors in both SL and TL regions is our next 
goal. Our model for the nucleon exploits many ingredients already introduced for the 
pion, adding suitable Ansatzes for the quark-nucleon amplitude in both the valence and 
nonvalence sectors. The possibility to address separately valence and nonvalence sectors 
of the hadron amplitude represents a distinct feature of our phenomenological model, and 
it is entailed by the model for the quark-photon vertex, we have adopted. In particular, 
within our approach one could ascribe the possible zero in the ratio G p E (Q 2 )fi p /G p M (Q 2 ) 
to the interplay between the contribution from the pair production process and the one 
from the elastic quark-photon vertex, shedding some light on the deep meaning of the 
puzzling experimental finding. 

The next step will be the extension of our calculation to the timelike region. 
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